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DEVELOPMENT OF A NONLINEAR UNSTEADY 


TRANSONIC FLOW THEORY 

By Stephen S. Stahara and John R. Spreiter* 

Nielsen Engineering & Research, Inc. 

SUMMARY 

A preliminary investigation was conducted to develop the basis of a 
nonlinear, unsteady, small-disturbance theory capable of predicting inviscid 
transonic flows about aerodynamic configurations undergoing both rigid body 
and elastic oscillations. The theory is based on the concept of dividing 
the flow into steady and unsteady components and then solving, by the 
method of local linearization, the coupled differential equation for the 
unsteady surface pressure distribution. 

The appropriate equations, valid at all frequencies, have been derived 
for two-dimensional flows at M = 1, and numerical results obtained for 
two classes of airfoils and two types of oscillatory motions. The results 
indicate that the theory smoothly and correctly converges to nonlinear 
quasi-steady theory as the reduced frequency of oscillation based on chord 
k 0 and to linear acoustic theory as k becomes large (k > 2) . Moreover, 
at low frequencies (k = 0.1) the theory shows quantitatively the large non- 
linear thickness effects induced on the unsteady flow by the steady motion 
and the complete inadequacy of linear theory in this frequency regime. 

INTRODUCTION 

The basic reason for developing an accurate unsteady transonic flow 
theory lies in the need to predict aerodynamic flutter and its associated 
dynamic instabilities, which are more likely to occur in the transonic 
regime than in any other. In view of emerging technological needs, however, 
such as those associated with the design of a transonic transport, the 
present state of unsteady transonic aerodynamics is clearly inadequate. 


♦Professor, Departments of Applied Mechanics and Aeronautics and Astro- 
nautics, Stanford University, Stanford, California. (Consultant at 
Nielsen Engineering & Research, Inc.). 



While some work has been done (see refs. 1-4) , such technically important 
problems as one-degree-of-freedom control surface flutter (control surface 
buzz) and the determination of stability derivatives at low reduced fre- 
quencies, to say nothing of the more general problem of predicting pressures 
on arbitrary oscillating surfaces at transonic speeds, remain essentially 
unsolved at the present time. Unless the frequency of oscillation is 
sufficiently high or the configuration of a special nature, agreement of 
existing theoretical results and data is not satisfactory (ref. 3) . The 
reason for these discrepancies stems from the linearization of the governing 
differential equations, which are fundamentally nonlinear in this speed 
regime. Consequently, it is clear that an adequate analysis of unsteady 
transonic flows must include at least some of the nonlinear effects present 
and one of the goals of this investigation is to identify and include the 
most important of these effects. 

Perhaps the primary reason for the relatively slow development of 
unsteady transonic aerodynamics is that, until recently, there existed a 
notable paucity of techniques for accurately predicting the steady tran- 
sonic flow about realistic aerodynamic shapes. Current research, however, 
has changed that situation so that the accurate theoretical treatment of 
realistic two- and three-dimensional steady transonic flows (refs. 5-10) 
is now possible. This is of major importance to the present technique 
which makes use of the nonlinear steady solution in determining the 
unsteady flow. 

Although the ultimate goal of this investigation is to evolve pre- 
dictive methods for calculating unsteady pressures, forces, and moments 
acting on general aerodynamic configurations throughout the transonic 
regime, the purpose of this preliminary study is to develop the basis of 
a nonlinear unsteady theory capable of accurately predicting transonic 
flows about airfoils undergoing both rigid-body and elastic oscillations. 


LIST OF SYMBOLS 


A 1 

equal to 

X^X, 

) 

a 2 

equal to 

(X 2 /2xW 1 - 4A 1 A 3 /A 2 2 


2 


c 



f 

F 

g 

X c 

k 

K 


Kl 


K- 


R 

R 


reference length: for two-dimensional flows, equal to chord 

length; for axisymmetric flows, equal to body length 

total pressure coefficient, equations (26) and (71) 
steady-state pressure coefficient, equation (25) 

unsteady pressure coefficient, equation (25) 

similarity form of steady pressure coefficient, equations (27) 
and (72) 

similarity form of unsteady pressure coefficient, equations (28) 
and (73) 

equal to dZ(x)/dx + ikz (x) 

function describing body surface, equation (11) 
equal to R (x) £dR (x) /dx + ikR(x)J 

zero-order modified .Bessel function of first kind 

reduced oscillatory frequency based on chord, ujc/U^ 

zero-order modified Bessel function of second kind 

similarity parameter, equal to - 1 /^M^t(y + 1)J 2/<3 

similarity parameter, equal to 2k/e (y + 1), equations (24) and 
(70) 

free stream Mach number 

dimensionless total body radius, normalized by c; equation (60) 

dimensionless affine function describing steady-state body radius, 
equation (60) 


R 

t 


dimensionless affine function describing amplitude of unsteady 
radial perturbation, equation (60) 

nondimensional time, normalized by c/O^ 


free stream— velocity 
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N i 


(x,y,z) nondimen sional , body-fixed Cartesian coordinate system with x 

axis directed rearward and aligned with longitudinal axis of 
body, y axis directed to the right facing forward, and z 
axis directed vertically upward; coordinates normalized by c 

x* location of steady-state sonic point, normalized by c 

Z dimensionless total wing or airfoil ordinates, normalized by c; 

equation (9) 

Zj dimensionless steady-state airfoil ordinates, normalized by c; 

equation (79) 

Z dimensionless affine function describing steady-state wing or 

airfoil ordinates, equation (9) 

dimensionless affine function describing amplitude of unsteady 
wing or airfoil ordinate perturbation, equation (9) 

y ratio of specific heats 

6 dimensionless amplitude of unsteady oscillation, normalized by c 

£ expansion parameter for steady flow, equations (16) and (63) 

expansion parameter for unsteady flow, equations (17) and (64) 

£ normalized similarity form of z ordinate, equation (18) 

X j equal to + K x 

X„ equal to ¥ + iK„ 

X 3 equal to - kK 2 /2 

5,^ dummy variables 

p normalized similarity form of radial ordinate, equation (65) 

t airfoil, wing, or body thickness ratio 

<j) dimensionless total perturbation velocity potential, equations (6) 

and (59) 

<J> dimensionless steady-state perturbation velocity potential, 

equations (6) and (59) 

$ amplitude of dimensionless unsteady perturbation velocity 

potential, equations (6) and (59) 
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$ similarity form of dimensionless steady-state perturbation 

velocity potential, equations (14) and (61) 

$ similarity form of dimensionless unsteady perturbation velocity 

potential, equations (15) and (62) 

A$ edge-correction potential to $ for subsonic flow, equation (58) 

to oscillatory frequency, rad/sec 


ANALYSIS 


Partial Differential Equations 

The concept that a major body of steady transonic flow problems can be 
predicted accurately within the framework of inviscid, small-disturbance 
theory without recourse to numerical solutions of the full nonlinear equa- 
tions has been established (refs. 11 and 12) for some time. There is good 
reason to believe that this premise also holds for the unsteady case, 
particularly for small-amplitude flutter and stability calculations. For 
unsteady motions, small-disturbance transonic theory yields the following 
differential equation for the dimensionless perturbation velocity potential 
<j> (ref. 3) 


(1 - Mi) 


XX 


+ <t> 


yy 


zz 


Mi, (y + 1) i 


X XX 


+ M 


°° T tt 


2M 2 4> . 
00 xt 


( 1 ) 


where is the free-stream Mach number, y is the ratio of specific 

heats, (x,y,z) are the spatial coordinates nondimensionalized by some 
reference length c, and t is the nondimen sional time normalized by 
c/U . 

Depending upon the nature of the unsteady motion, equation (1) can 
simplify into various forms. For very fast oscillations, the following 
linear, two-dimensional wave equation holds. 


yy 


zz 


= M 2 


K tt 


(2) 


while for slightly slower unsteady motions, the more complicated linear 
equation is valid. 

(1 - M 2 )c(> + <(> + <p = M 2 <j> ,, + 2M 2 cj> . 

°° T xx yy T zz =° T tt °° T xt 
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For slowly unsteady motion, the following nonlinear equation applies, 

(1 - M 2 )(f> + <J> + $ = M 2 (y + 1) <J> 4> + 2M 2 d> . (4) 

00 xx yy zz °° Y x xx °° y xt ' ' 

while for very slow (quasi-steady) motions, the same nonlinear equation 
that holds for steady flow applies. 

(1 - M 2 ) <j> + <j> + <j> = M 2 (y + 1 ) 4> d> (5) 

Solutions to the linear two-dimensional wave equation (eq. (2)) are 
straightforward and not of primary concern here. Miles (ref. 13) has 
treated the linear equation (3) by analytical means for supersonic flows, 
while Landahl (ref. 1) has applied the same equation at M = 1 to certain 
classes of wings and low aspect ratio wing-body combinations. While many 
important applications remain to be worked out for this equation, it is 
clear from past comparisons with experimental data that some of the non- 
linear features characteristic to problems in this speed regime, which have 
been omitted by this equation, must be included. 

Very little work has been done on solving the more difficult nonlinear 
equation (4). In many practical flutter and stability problems, the reduced 
frequencies are very low and in investigations of dynamic stability, the 
values of aerodynamic stability derivatives as the reduced frequency k ■* 0 
are required. Furthermore, it can be shown by the method of matched asymp- 
totic expansions that equation (4) provides the proper solution for a 
relatively wide range of low frequency unsteady motions and that the solu- 
tions of equation (4) converge uniformily as k -*• 0 to the nonlinear 
quasi-steady results provided by equation (5). Thus, solutions to equa- 
tion (4) contain information about both low frequency and quasi-steady 
nonlinear motions. Consequently, while in this study we will seek solutions 
to the more general nonlinear unsteady small-disturbance equation ( 1 ) , we 
note that for low frequency motions the simpler nonlinear equation (4) 
could be legitimately used. 

Finally, with regard to steady-state solutions of the nonlinear equa- 
tion (5), several comments are appropriate. First, it has been adequately 
demonstrated (refs. 11 and 12 ) that solutions of this equation are capable 
of providing good representations of actual steady transonic flows about a 
wide variety of shapes. Applications to thin airfoils (refs. 14, 15, and 
16) finite-span wings (ref. 17), slender bodies of revolution (refs. 15 f 
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18, and 19), slender nonaxisymmetric bodies (ref. 8), and even certain 
classes of wing-body combinations (refs. 8, 9, and 10) have confirmed this. 
Consequently, the nonlinear effects associated with a major body of tran- 
sonic flows can be adequately accounted for by inclusion of the 't’x'^xx 
term in equation (5) . Furthermore, at this time both exact and approxi- 
mate methods for solving this equation are available as indicated above 
(ref s . 5 , 6 , and 7) . 


General Oscillatory Motions 

Fundamental to any study of unsteady aerodynamics is the analysis of 
oscillatory motion. For general oscillatory motions, it is convenient to 
expand the solution into a steady and unsteady component. Thus, we may set 

<f>(x,y,z,t) = + R.P. <f> (x,y,z)e lkt ^ (6) 

where k is the reduced frequency of the oscillatory motion, R.P. signi- 
fies the real part of the complex quantity, <J> is the steady perturbation 
potential, and <J> is the amplitude of the oscillatory potential. Inserting 
equation (6) into the full equation (1) , the following coupled partial 
differential equation is obtained for <j> . 


(1 - M*)* 


xx 


yy 


zz 


M (y + 1) (4> 4> + 4> , 4> 

“> x xx ^ i xx 

x 


x / 

XX 


M 2 k 2 d> + 2iM z kd> 

on T co T ^ 


(7) 


This equation has variable coefficients, is both nonlinear and of mixed 
elliptic-hyperbolic type and is, in general, more difficult to solve than 
the corresponding one for the steady component (eq. (5)). Since, in 
flutter and stability analysis, the stability of only small-amplitude 
oscillations is normally investigated, it is appropriate to assume that 
the oscillatory flow is a small perturbation on the nonlinear steady flow. 
This is logical inasmuch as oscillations or other unsteady deformations of 
the surfaces of an aircraft should involve only minor deviations from their 
steady-state positions. With this simplification, the governing equation 
for <j) becomes 
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(1 - M*)' 


XX 


yy 


zz 


= M 2 (y. + !)(<(>, 4) 


x 


XX 


XX 


+ 2iM kd> - 

oo Y x 


M 2 k 2 ( 


(8) 


which, although being linear remains quite formidable to solve because of 
the variable coefficients and mixed elliptic-hyperbolic type. It is likely, 
however, that solutions of this equation will provide all of the infor- 
mation necessary for an accurate transonic stability analysis. 

With the above approach of splitting the perturbation potential, the 
surface boundary condition can be treated in a very general way, accounting 
for both rigid body and elastic oscillations. For example, for thin wings 
where Z(x,y,t) represents the normalized ordinates of the upper surface, 
we can set 


Z (x,y , t) = tZ (x,y) + R.P. ^6Z (x,y) e lkt ^> (9) 

where t and 6 represent, respectively, the wing thickness ratio and the 
normalized amplitude of the unsteady oscillation. In this representation 
rigid body plunging oscillations can be given by Z(x,y) = 1, pitching 
oscillations about a line at x - x Q by Z(x,y) = (x - x ) , elastic longi- 
tudinal oscillations by Z(x,y) = f (y) sin (mrx) , etc. Solutions for more 
general unsteady motions such as gust reponse, buffeting, etc., can be 
obtained by superposition of such oscillatory solutions. 

The surface boundary condition for the oscillatory component is found 
from the expression 

3F 

— + grad (x + ij))*grad F = 0 on F(x,y,z,t) = 0 (10) 


where 


F(x,y,z,t) = z - Z(x,y,t) 

with <p given by equation (6) and Z(x,y,t) by equation 

+ ilcZ(x,y)j 

At infinity, the unsteady velocity components must vanish 
fashion. 


(x,y,0±) = ±6 


3 Z (x , 
3x 


(ID 

(9) . Thus , 

( 12 ) 

in an appropriate 


■cf* 
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In transonic and supersonic flows, it is also necessary, in general, 
to provide appropriate relations for the discontinuous changes in velocity 
that occur at shock surfaces. However, in the case of sonic and near-sonic 
flows which are considered in this report, for large classes of thin air- 
foils and wings with convex contours at zero or small angles of attack, the 
shock waves present originate at the trailing edge. For small amplitude 
oscillations, they remain substantially fixed at that location and do not 
require further attention. However, we note that in general even for these 
classes of shapes, at Mach numbers in the lower transonic regime when shock 
waves appear on the upper and/or lower surface and at Mach numbers in the 
upper transonic regime when bow shocks appear, these shock relations must 
be properly accounted for. 

Finally, the pressure coefficient for the thin wing case is given by 


C P = 


R.P.{^ X + ik^ 


ikt 


(13) 


while for slender bodies additional terms, as usual, appear (see Similarity 
Equations for Axisymmetric Flows) . 


Similarity Equations for Two-Dimensional Flows 

By using the method of matched asymptotic expansions, it can be shown 
that for two-dimensional flow the potentials (<P lt <J> ) can be expressed in 
the similarity form 

cf> i (x,z) = e?(x,C) + 0 (e 2 ) (14) 


<t> (x,z) = e 1 $(x,C) + 0 ( e 2 ) 


(15) 


where 


Lm^(y + 1) j 

[tm 2 (y + 1)] 1/3 


(16) 


(17) 
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C = [m*(y + 1) T j 17,3 z 


and the similarity potentials ($, $) satisfy the following differential 
equations and boundary conditions 


-Fk + $ I $ + 

L 1 x J xx 


V ‘ 0 


V x ' 01) , t wp. 


-[k + $ 1 $ + i „ = [$ + iK 1 

Li xJxx SC Lxx 2 J 


$ x - y V 


- 

:, 0 ±) = ± 

. 


dZ (x) 


+ ikZ(x) 


where 


M* - 1 


1 [m^t (Y + 1)1 


K = 2 £ = 2k ” 

2 y + 1 e t (y + 1 ) 


The physical pressure coefficient can be represented by 


C = C + R.P.fc e 

P Pi V p 


with C defined by equation (13) , or in similarity form by 
P 


C = ec + e,R.P .f 

P P 1 V 


C e ikt ■) 
P J 
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where and are the similarity forms of the steady and unsteady 

pressure coefficients and are given by 


C = - 2 $ 

P x 


(27) 


C = -2 f $ + ik$"\ 

p V x J 


(28) 


Thus, the differential equation and boundary condition for the unsteady 
component can be expressed in the compact form 


A $ + A 4> + A , $ — 4> 

1 XX 2 x 3 C £ 


(29) 


4>^(x,0±) = ± f ( x ) 


(30) 


where 


K, + $ 

1 x 


4> + iK, 

XX 2 


(31) 

(32) 


. I K 

2 K 2 


(33) 


f (x) = dZ d ( x x) + ikZ(x) 


(34) 


so that, depending upon the value of Aj, the following set of equations 
must be solved s 


(1) Sonic Equation 


A 2 4> x + A 3 * = 


(x = x*, Aj = 0) 
(R.P. {A 2 } > 0) 


(35) 


(2) Supersonic Equation 


(x > x*, A 3 > 0) 


A i $ xx + A 2 $ v + A 3 4> - $ 


x 




(36) 


(3) Subsonic Equation 


(x < x*, A < 0) 


^ X l )$ xx - X 2 $ x - + = 0 


(37) 
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where x* denotes the location of the steady flow sonic point, i.e. , 
where A = 0 . 

The occurence of the three fundamentally different equations for $ 
is necessary and is a result of the unsteady flow disturbances riding 
on the basic characteristics of the steady flow which contains subsonic 
(Aj < 0), sonic (Aj = 0), and supersonic (A i > 0) regions. The regions 
in which the various solutions will apply can be seen more clearly in the 
figure below, which represents schematically the typical steady surface 
pressure distribution on a thin, convex airfoil at = 1. 



The solution for the sonic region (2) must merge continuously with that for 
the subsonic region (1) ahead of it and the supersonic region (3) behind it. 

With regard to the A 3 4> term which appears in the preceeding equa- 
tions, if we were to restrict the analysis to slowly unsteady nonlinear 
motions, then this term, which results from the M^<|> t term in equation (1) , 
could be legitimately omitted (cf. eqs. (1) and (4)). However, it appears 
that no undue complications arise from retaining this term. Consequently, 
the analysis and theory presented in this report are applicable, without 
restriction, to all frequencies. 

Local Linearization Method 

The local linearization method , originated by Spreiter and Alksne 
(refs. 14, 17, and 18), was developed as an approximate method for solving 
the nonlinear transonic steady flow equation (5) , although the basic ideas 
of the technique are apparently applicable to a large number of both linear 
and nonlinear problems. It has been applied successfully in the past to 
transonic flows past a wide variety of thin airfoils (refs. 14, 15, and 16), 
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slender bodies of revolution (.refs. 15, 18, and 19) and wings of finite 
span (ref. 17), and more recently to even more complex shapes, including 
nonaxi symmetric bodies and wing-body combinations (refs. 8, 9, and 10). 

With regard to the present problem, the procedure for determining the 
local linearization solution for the unsteady component $ is to replace 
temporarily the variable coefficients (X x , X ) by constants, solve the 
simplified equations (35), (36), and (37) for 4>(x,£), calculate the un- 
steady surface acceleration d 2 $ (x , 0 ) /dx 2 , then replace the constants 
( X x , X ) by the functions they originally represented, and finally integrate 
the resultant second-order ordinary differential equation to obtain both 
$^(x,0) and $(x,0) as are required (see eq. (13)) to determine the unsteady 
surface pressure. 

Solutions of equations (35) , (36) , and (37) with (X 1 , X ) as con- 

stants have been determined and are given for the sonic case (X x 2 0) by 





sonic 



f (C)e 


M 2 

' (x ?) 4 (x-C) 


V x - s 




(38) 


for the supersonic case (X > 0) by 


($ (x, O') 


super yr 
sonic » i o 


X 

kS 


- Vm 


f (£)e 


- Aj (x-C) 


A 2 (x - £) z - X^j 


d£ (39) 


and for the subsonic case (Xj < 0) by 

1 

= ------- f 

TV V(-X, ) J 


u) 

sub 


sonic 

A ! 

and 


fU)e 


- Aj (x- 5 ) 


K 


^ - A 2 -^(x - o 2 + (-X t ) c 2 j 


d? 

(40) 


A i 2 X . 


A 2 2X 


.1 

■ V 1 ' — 


(41) 


(42) 
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and where I and K Q are zero-order modified Bessel functions of the 
first and second kind, respectively. As x x* (i.e., as X x -*■ 0) , it 
can be shown that the supersonic and subsonic solutions (and their deri- 
vatives) given by equations (39) and (40) uniformily converge to the sonic 
result given by equation (38) as they must. 

Application of the local linearization method, together with some 
manipulation to obtain forms convenient for numerical integration, provides 
the following second-order ordinary differential equations for the unsteady 
surface acceleration d 2 0 (x , 0) /dx 2 . 


/ d 2 i> (x, 0) \ 

V dx 2 ) 


-1 


- f (0)e 


sonic 


Y7T 2 L V*~ 


(k + xj} 


f ' ( 0 ) e 


*7' 




+ 2f"(x) e 


r- A 3 

V* - x^ * {f"(x) - f"(?)} e 


- x^x-5) 


/• 1 de '/ 


V x - s 


d? 


(43) 


f d 2 $ (x , 0) \ 
\ dx 2 J 


super 

sonic 


-1 

vt; 


- Ax 

f (0)e -A 


T o [ - xf- - *■ [*»*]} 


- Ax 

+ f ' (0)e -I 


M + b 


(C)e 


- A, (x-C) 


A 2 (x - ?) 


H 


(44) 


/ d 2 $ (x , 0 ) \ 
\ dx 2 ) 


-1 


. ttV-X 

sub » l 


sonic 


- f (1) e 


-Ax 

f (0) e 1 • (— Aj ) • 


A, (1-x) 


{ k o[' a = x ]-S 7 k .[- r ^]} 


• <'V f K o [ - A, U - *>] + ^ K,[- <1 - *>] J 

-Ax r A 

+ f ' (0)e • K o - A 2 x - f ' (1) e 
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(1-x) 

(Continued on next page) 


K o [-A 2 (1 -x) 


fl 



+ / {£ " 


(?) 


- f"(x)} 


- A, (x-5) 


K 


A 2 (x - 5) j 


d? 


1 

+ J £f"(?) 

X 



A x (?-x) 
e 




d? 



where now 
i.e. 





X ) and, consequently, (A 



, A ) are functions of 


x 


(45) 
only , 


A j = 4> x (x,0) + K x (46) 

A 2 = ^ xx (x , 0) + iK 2 (47) 

In order to integrate equations (43) , (44) , and (45) , two boundary 
conditions are required. These are provided, in direct analogy to the 
local linearization solutions for the steady case, by requiring that at 
the steady-state sonic point x* 
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-1 


$(x*,0) = ($(x*,0)) 


sonic 


V 


== f 
"V 4 


f (g)e 


2 


V x * - ? 


d6 


(48) 


_ Z' d$ (x* , 0 ) 
“V dx y 


d$ (x* , 0 ) _ f d$ (x* , 0 ) ' 
dx 


-1 


sonic 




I 


f (0)e 


V 


x’ 


v** 


f' me 


V 


(x*-5) 


V x * - ? 


dC 


(49) 


where X 2 * = X 2 (x*) = $ xx (x*,0) + iK 2 . 

For numerical evaluation, more convenient forms of equations (48) and 
(49) are 


$ (x* , 0 ) 


-1 

V^v* 


Vr - 


2f (x*) 


/ 






(x*-C) 


V"x* - e 



(50) 


d$ (x* , 0 ) 
dx 


-1 


f (0)e 


*2* 


Yx* 



o 


e 






y** - s 


(x*-C) 



(51) 


Although in their general form, equations (43) , (44) , and (45) appear 
formidable, for certain oscillatory motions, they simplify considerably. 

For example, for plunging oscillations where Z (x) =1 so that f (x) = ik, 
we have 
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Edge Correction for Subsonic Potential 

Because of the absence of upstream influence, the sonic and supersonic 
solutions given above are complete, and no further boundary conditions for 
the wake or for points ahead of the airfoil are required. For the subsonic 
solution, however, these additional conditions must be taken into account. 

For purely subsonic (A l < 0 everywhere) unsteady lifting flows, these 
requirements on 9 are shown in the figure below 

z 

x , 0) = 0 

x 

A j < 0 

where the condition that $(x,0) = 0 ahead of the airfoil insures the 
physical requirement that no discontinuity occurs in 9 when crossing the 
z = 0 line, and the condition that 9 (x,0) + ik5(x,0) =0 in the wake 
insures the continuity of pressure along that line. For a supersonic 
trailing edge, which is the situation for the - 1 flows considered 

herein, the wake condition can be dropped. To satisfy the condition that 

<I>(x,0) = 0 ( - <» < x < 0) (55) 
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an additional potential, 
tion (40) such that A4> 
cancels the effect of $ 


A4>, must be added to $ su bsonic 9^ ven ky equa- 
satisfies the subsonic differential equation (37) 

subsonic ahead of the airfoil, 


A$(x,0) 


^subsonic 


(x,0) 


( - 00 < x < 0) 


(56) 


and does not alter the boundary condition satisfied by $ . on the 

J subsonrc 

airfoil surface. This can be accomplished by requiring that (see ref. 1) 

A$ z (x,0) =0 (0 < x < “ ) ( 57 ) 

With these conditions, the solution for A 5 has been found and on the 
x-axis is given by 


A<f> (x, 0) 


1 

* 2 V7-*7> 


/ 


- A, (x-5) 

f (?) e 


I 



e 


Aj (x+C,) 
x + ?! 


X 


K 


o 





d? 2 dC 


(58) 


It can be shown that in the vicinity of the steady sonic point (i.e., 
as x -*■ x* , X j -*■ 0) , the contribution of A4> (and its derivatives) to the 
unsteady pressure distribution uniformity go to zero. Consequently, the 
major contribution of A4> occurs in the vicinity of the leading edge, so 
that A4> can be regarded as essentially an edge-correction. We note that 
in the numerical results presented in this report, A$ has not been 
included . 


Similarity Equations for Axisymmetric Flows 

An analysis parallel to that used for the two-dimensional case has 
been carried out for axisymmetric flows in order to establish in that 
situation the governing equations and boundary conditions for the unsteady 
component. By splitting the potential into a steady and unsteady component, 

<j> (x,r ,t) = <(> j (x,r) + R. P . ^ (x ,r ) e ikt ^ (59) 
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and the normalized body ordinates in a similar fashion. 


l(x,t) = tR(x) + R.P. ^ 6R(x)e ikt ^ 


and then using the method of matched asymptotic expansions, the potentials 
(<!> j , 4>) can be given in the following similarity form 

4>! (x,r) = £^(x,p) + 0(e 2 ) (6i; 

‘l’ (x,r) = e x $(x,p) + 0 ( £ 2 ) (62) 

where 


£ x = X 6 


p = TM V(Y + 1 ) r 


and the similarity potentials ($, 3>) satisfy the following differential 
equations and boundary conditions 


-[ K > + *J 


$ + $ = o 

xx p p p p 


lira (Pfc) = R (x) 

p+0 p ax 


-Tk. + 0 $ + $ = [$ + iK 1 

L 1 x J xx pp pp L xx 2 J 


$ - 9- K,<i> 

x 2 2 


lim f pi "V = R(x) -A -—- - + ikR(x) 

p-*0 V 0J 1 dx 


where K x is given by equation (23) and 


K = 2 k = 2 jc_ 

2 Y+le Y+1 t 2 


The physical pressure coefficient in this case is given by 


c p * + e 


. s - p -(v ikt ) 
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where (e, e 1 ) are given by equations (63) and (64) and the similarity 

pressure coefficients (C , C } are defined as 

P P 

C = - 2 T - F 2 (72) 

p x p 

C = - 2 r i + ikO + $ 5 I (73) 

p L x p p J 

Thus, the differential equation and surface boundary condition for 
the unsteady component $ can be expressed in the form 


X i $ + X $ + 

1 XX 2 X 


X $ 

3 


$ + i $ 

PP P P 


(74) 


lim (p$ ) = g(x) (75) 

P+0 p 

with Xj , X 2 , and X 3 defined as before by equations (31), (32), and (33) 
and g (x) given by equation (69). 

The solutions to equation (74) for constant (X , X 2 ) needed as the 
first step in the local linearization method have also been determined. 

They are 


x 

( i(x ' p) ) sonic = - 2 - f 


g ( £ ) exp 


(x - E.) 
x - ? 


X 2 P Z 

4(x - S) . 


d? 


(76) 


(A, = 0) 


( i(x ' p) ) super = - / 
sonic ~ 




g (?) x 




exp -xr — (x - ?) 


cosh 


4X^3 


2X . 


A 2 


V(x -?) 2 - Xji 


V(x - ?) 2 - Xj 


- d? 


(Xj > 0) 


(77) 
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O' (x, p)) 


sub 

sonic 


if 


g (?) x 


exp 


(~2X 1 ) 


x - C — v / 1 - 


4X^3 
*2 


V(x - C) 2 + ( 


■ x ‘ )p! } 


Vu - O 2 + (-X ! ) 


d5 


(X 1 < 0) 


(78) 


RESULTS AND DISCUSSION 

In order to illustrate the general behavior of the unsteady results 
predicted by the present theory, the two-dimensional unsteady local lineari- 
zation equations (43) , (44) , and (45) have been programmed and results 

obtained for two types of oscillatory motions for two different classes of 
airfoils at =1. A member of the first class of airfoils (Guderley 

airfoils) is shown schematically in figure 1 together with its steady-state 
similarity pressure distribution C for = 1 . These airfoils, whose 

normalized ordinates Z 1 are given by (ref. 14), 

Z j (x) = t(x) 3/2 (1 - x) (79) 


are characterized by a constant steady-state surface pressure gradient 
8Cp^/3x in the transonic Mach freeze range, i.e., M m - 1. Their similarity 
steady surface pressure distribution at =1 is given by 


C 

P 


(Y + 1) 


1/3 


TP 




(80) 


The second class of airfoils, shown in figure 2, is composed of parabolic- 
arcs with normalized ordinates given by 


Z j (x) = 2x (x - x 2 ) 


(81) 
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and with similarity surface pressure distribution at M = 1 represented 
accurately by the local linearization result (ref. 14) 


= - 2 


i r ( ln [ 4x ] " 8x + 8x2 + I) 


1 V 3 


(82) 


The two types of oscillatory motions considered are vertical translation 
(Z = 1) and pitching about the nose (Z = x) . 

Figure 3 exhibits the solutions for the normalized magnitude and phase 
of the unsteady surface pressure distributions on the upper surface of a 
6-percent thick Guderley airfoil oscillating in vertical translation 
(plunging) at various reduced frequencies k at = 1. Note that in 

order to provide a form more convenient for comparison, in all of the 
results presented herein for plunging oscillations the magnitude of the 
unsteady pressure distribution has been normalized by 5k, rather than the 
usual ' 6 since in this case $ contains a linear factor in k (see 
eqs. (52), (53), and (54)). The results of figure 3 illustrate the pro- 
gression, for a fixed airfoil geometry as k goes from large values to 
zero, of the results of the present theory as they converge smoothly to 
the nonlinear quasi-steady result. In particular, we note that for small 
values of k, i.e., k - 0.1, the present results are close to those of 
quasi-steady theory, implying that in the low frequency regime, the non- 
linear thickness effects of the steady flow are of primary influence on 
the unsteady flow and that, for very low frequencies, quasi-steady theory 
can provide a reasonable approximation. 

While figure 3 demonstrates the ability of the present theory to 
converge uniformily to the nonlinear quasi-steady result as k -*■ 0 , 
figure 4 shows the corresponding ability, as k becomes large, of the 
present theory to converge to linear acoustic theory results. For low k 
(i.e., k = 0.1), figure 4 clearly indicates that linear theory is virtually 
useless (the phase angle for k = 0.1 is off the scale of the plot) , while 
at larger values (k = 2) , the results of the present nonlinear theory and 
linear acoustic theory are converging. At k = 5, the results are 
essentially identical. 

Figures 5 and 6 show the analogous results of figures 3 and 4 for a 
6-percent thick parabolic-arc airfoil. Figure 5 demonstrates that, as 
was the case for the Guderley airfoil, the present theory moves smoothly 


22 



toward the quasi-steady result as k becomes small, while figure 6 
illustrates the corresponding convergence to linear theory as k becomes 
large. In particular, we note in comparing figures 4 and 6, that as k 
becomes large, the results for both the Guderley airfoil and the parabolic- 
arc airfoil uniformily approach each other (as well as linear theory) , 
which agrees with the concept that for large k the unsteady motion 
becomes independent of the profile geometry. 


In figure 7(a) and (b) , the normalized magnitude (|c | / <5 ) and phase 

1 ? 

of the unsteady pressure distribution are presented for a 6-percent thick 
Guderley airfoil oscillating in pitch about the nose at M m = 1. For this 
unsteady motion, we note again the smooth variation of the present theory 
to the quasi-steady result as k -*■ 0 . Figure 8(a) and (b) compares these 
same results with linear theory. This figure demonstrates that for pitching 
oscillations, as was the case for plunging motion, linear theory is in 
serious error for Jc less than 1. For k = 2, the magnitude of pre- 

dicted by linear theory is in good agreement with the present nonlinear 
results while the phase is only in fair agreement. At k = 5, both magni- 
tude and phase are essentially identical to the present theory. Figures 
9(a) and (b) and 10(a) and (b) show the corresponding results for a 6- 
percent thick parabolic-arc airfoil. These demonstrate virtually the same 
behavior as those shown for the Guderley airfoil. 


Since all the results presented thus far have been for a fixed thick- 
ness ratio, in order to demonstrate explicitly the importance of nonlinear 
thickness effects we have prepared figure 11 which shows, for a fixed 

reduced frequency (k = 0.1), the normalized magnitude (|C | /6k) and phase 

P 

of the unsteady surface pressure due to plunging oscillations of various 
thickness-ratio parabolic-arc airfoils. We note the enormous change in 
both magnitude and phase as x increases only slightly from zero (for 
which the present theory and linear acoustic theory are the same) to one- 
percent thickness ratio (x = 0.01). These results demonstrate quantita- 
tively both the importance of thickness effects and the inadequancy of 
linear theory at low frequencies, even for exceptionally thin airfoils. 
Figure 12 shows similar results for the higher reduced frequency k = 1.0. 
In this case, the differences in the magnitude of the unsteady pressures 
for the various thickness ratios are not nearly as great as for k = 0.1 
and also are largely restricted to the aft portion of the airfoil. The 
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phase angles, which are generally more sensitive to change remain somewhat 
separated for each thickness ratio. Nevertheless, the trend of all the 
curves to move toward the zero thickness (linear) result is clear. This 
trend is notably demonstrated in figure 13 which shows the analogous 
results for k = 10. In this case, there is essentially no difference in 
the magnitude or phase angle for the various thickness ratios from the 
linear (t = 0.0) result, except near the tail, where presumably the flow 
finally has had enough time to react to some of the nonlinear effects of 
thickness in the steady flow. 


CONCLUDING REMARKS 


Theoretical analysis has been carried out to develop the basis of a 
nonlinear, unsteady, small-disturbance theory valid for all frequencies 
which is capable of predicting transonic inviscid flows about thin airfoils 
undergoing both rigid body and elastic oscillations at M - 1. The theory 
employs the concept of dividing the flow into steady and unsteady components 
and then solves, by the method of local linearization, the coupled unsteady 
equation for the surface pressure distribution. 

Calculated results for two classes of airfoils and two types of 
oscillatory motions indicate smooth convergence to the quasi-steady results 
as the reduced frequency based on chord k ->■ 0 and to linear acoustic 
results as k becomes large (k > 2). Moreover, at low frequencies 
(k = 0.1) the theory demonstrates quantitatively both the large nonlinear 
thickness effects induced on the unsteady flow by the steady motion and the 
subsequent inadequacy of linear theory in this frequency regime. 

Application to more complex aerodynamic configurations is possible, 
as shown by the derivation and partial solution of the corresponding 
results for unsteady axisymmetric flows. At this stage of development, 
several details of the theory remain to be worked out. For example, the 
proper unsteady shock relations need to be determined and applied, both 
to surface and bow shocks, and several questions concerning the error 
incurred in the unsteady surface pressure by using the local linearization 
approximation need to be investigated. Experimental verification of the 
theory is essential at this stage and while we are unaware of any data 
suitable for comparison with the preliminary results presented herein, 
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future work should be directed toward those configurations for which 
experimental data exists. Finally, the results obtained in this prelim- 
inary investigation strongly suggest that further development and analysis, 
to extend and refine this promising technique, be carried out as soon as 
possible in order to provide a practical and accurate aerodynamic tool for 
analyzing unsteady transonic flutter and stability problems. 


Nielsen Engineering & Research, Inc. 
Mountain View, California 
November, 1972 
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Figure 1 . -Similarity form of steady surface 
pressure distribution for a family of 
affinely-related Guderley airfoils at 
M =1 according to the method 
of local linearization. 



Figure 2 . -Similarity form of steady surface 
pressure distribution for a family of 
affinely-related parabolic-arc airfoils 
at M ro = 1 according to the method of 
local linearization. 
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6 . -Comparisons of local linear 
steady pressure distributions g 
in Figure 5 with linear theory 






Figure 7 (a) . -Magnitude of unsteady pressure distributions 
on the upper surface of a 6-percent thick Guderley 
airfoil oscillating in pitch about the nose at 
various reduced frequencies k at M = 1. 




Figure 7(b). -Phase of unsteady pressure distributions 
on the upper surface of a 6-percent thick Guderley 
airfoil oscillating in pitch about the nose at 
various reduced frequencies 15 at M = 1. 
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Figure 8 (a) .-Comparisons of the magnitude of the 
local linearization unsteady pressure 
distributions given in Figure 7 (a) with 
linear theory. 
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Figure 8 (b) .-Comparisons of the phase of the local 
linearization unsteady pressure distributions 
given in Figure 7(b) with linear theory. 




Figure 9 (a) . -Magnitude of unsteady pressure distributions 
on the upper surface of a 6-percent thick parabolic-arc 
airfoil oscillating in pitch about the nose at various 
reduced frequencies Jc at M = 1 . 





Figure 9(b). -Phase of unsteady pressure distributions on 
the upper surface of a 6 -percent thick parabolic-arc 
airfoil oscillating in pitch about the nose at 
various reduced frequencies k at M =1. 







Figure 10 (a) .-Comparisons of the magnitude of the local 
linearization unsteady pressure distributions 
given in Figure 9(a) with linear theory. 
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Figure 10 (b) . -Comparisons of the phase of the local 
linearization unsteady pressure distributions 
given in Figure 9(b) with linear theory. 
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Figure 11 . -Magnitude and phase of unsteady pressure 
distributions on the upper surface of parabolic-arc 
airfoils having various thickness ratios oscillating 
in vertical translation at the reduced frequency 
k = 0.1 with M =1. 
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Figure 13 .-Magnitude and phase of unsteady pressure 
distributions on the upper surface of parabolic-arc 
airfoils having various thickness ratios oscillating 
in vertical translation at the reduced frequency 

k = 10. with M =1. 
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